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Abstract. H 1504+65 is an extremely hot hydrogen-deficient white dwarf with an effective temperature close to 200 000 K. We 
present new FUV and soft X-ray spectra obtained with FUSE and Chandra, which confirm that H 1504+65 has an atmosphere 
primarily composed of carbon and oxygen. The Chandra LETG spectrum (60 - 160 A) shows a wealth of photospheric absorp- 
tion lines from highly ionized oxygen, neon, and - for the first time identified in this star - magnesium and suggests relatively 
high Ne and Mg abundances. This corroborates an earlier suggestion that H 1504+65 represents a naked C/O stellar core or 
even the C/O envelope of an O-Ne-Mg white dwarf. 
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1. Introduction 

H 1504+65 is a faint blue star that has been identified as the 
counterpart of a bright soft X-ray source (Nousek etal. 1986) 
discovered by an early X-ray survey (Nugent etal. 1983). 
Spectroscopically, the star is a member of the PG1159 class 
which comprises hot hydrogen-deficient (pre-) white dwarf 
stars (r eff =75 000K-180000K, logg=5.5-8 [cgs]; Werner 
2001). The PG1159 stars are probably the outcome of a late 
helium-shell flash, a phenomenon that drives the currently 
observed fast evolutionary rates of three well-known objects 
(FG Sge, Sakurai's object, V605 Aql). A late helium-shell flash 
may occur in a post-AGB star or a white dwarf. Flash induced 
envelope mixing generates a H-deficient surface layer. The sur- 
face chemistry then essentially reflects that of the region be- 
tween the H-and He-burning shells of the precursor AGB star. 
The He-shell flash transforms the star back to an AGB star and 
the subsequent, second post-AGB evolution explains the exis- 
tence of Wolf-Rayet central stars of planetary nebulae and their 
successors, the PG1 159 stars. 

Within the PG1159 group H 1504+65 is an extraordinary 
object, as it has been shown that it is not only hydrogen- 
deficient but also helium-deficient. From optical spectra it was 
concluded that the atmosphere is primarily composed of car- 
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bon and oxygen, by equal amounts (Werner 1991, W91). Strong 
neon lines were detected in soft X-ray spectra taken with the 
EUVE satellite and in an optical-UV Keck spectrum and an 
abundance of Ne=2-5% (all abundances in this paper are given 
as mass fractions) was derived (Werner & Wolff 1999, WW99). 
The origin of this exotic surface chemistry is unclear. 

Here we present new results of observations performed 
with FUSE (Far Ultraviolet Spectroscopic Explorer) and the 
Chandra X-ray Observatory, whose spectroscopic resolution is 
more than an order of magnitude better than that of previous 
FUV and EUV missions (HUT and EUVE). We will show that 
the Chandra spectrum contains a wealth of photospheric ab- 
sorption lines from highly ionized metals. Ionization balances 
of O, Ne, and Mg lines provide new constraints on the effec- 
tive temperature and, for the first time, allow an estimation of 
the Mg abundance. We also analyze the spectra with respect to 
other metals, namely Al, Na, and the iron group. 

In the following sections we will first describe the model at- 
mosphere calculations. Then we will present in detail the com- 
pilation of atomic data and design of the model atoms. After 
a coarse characterization of the FUSE and Chandra spectra we 
turn to the detailed comparison between synthetic and observed 
spectra. Finally we discuss the results and their implications on 
the evolutionary history of H 1504+65. 

2. Model atmosphere calculations 

Line blanketed NLTE model atmospheres were computed us- 
ing our PR02 code (Werner et al. 2003a). The models are in 
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Table 1. Summary of the model atoms used in the model atmo- 
sphere and line formation calculations. The numbers in brack- 
ets give the individual line numbers summed into superlines for 
the heavy metal ions. 
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Fig. 1. Chandra observation (thin line) compared to two model 
spectra (thick lines). From one of the models we have omit- 
ted the bound-free cross-section of the level Mg v 2p 5 3 P° in 
order to demonstrate that a strong autoionization feature in 
the cross-section (dashed line) can result in a line-like absorp- 
tion feature in the calculated spectrum (marked by the arrow, 
r eff =175 000K). 

hydrostatic and radiative equilibrium. The atomic models ac- 
count for the most abundant elements in H 1504+65, carbon, 
oxygen and neon. Helium is also included in order to be able 
to derive an upper limit for the He abundance from the He n 
Balmer lines in the FUSE spectra. It is important to include 
Stark broadened profiles for the strongest spectral lines (from 
C iv, O vi and Ne vn) even at this stage of calculations, because 
they do affect the atmospheric structure by blanketing and cool- 
ing effects (Werner 1996). 

Model atmospheres with numerous parameter sets for T e ff, 
log g, and He/C/O/Ne element mixtures were computed to find 
the best fitting model for H 1504+65. Based on these model 
structures (keeping them fixed), we performed NLTE line for- 
mation iterations for additional chemical elements (Mg, Al, 
Na, Fe, and Ni) in order to guide our search for lines of these 
species in the observed spectra. Test calculations have shown 
that neglecting the back-reaction of these species onto the at- 
mospheric structure and the emergent spectrum is unimportant 
for our analysis, except for the iron group elements, which 
can affect the overall flux distribution when included into the 
models self-consistently. Hence, we have also computed ex- 
ploratory models accounting for all elements with Z=20-28 
(Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni), using a generic model atom 
comprising these species. Table[0summarizes all model atoms 
employed. For the final spectrum synthesis fine structure split- 
ting of atomic levels of light metals must be considered; level 
populations were taken from the models with appropriate sta- 
tistical weights. For some multiplets of neon and magnesium, 
treatment of the fine structure splitting is not entirely possible 
because excitation energies are not known for all sublevels. As 
a consequence, uncertainties of the order 0. 1 A in line positions 
within a multiplet remain. 

In the following section we describe in detail the design of 
our model atoms used for the spectral analysis. This is neces- 
sary and useful because our Chandra spectrum of H 1504+65 
represents the first high-resolution spectral observation of a 
purely photospheric stellar spectrum in the soft X-ray region. 
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Table 2. Photospheric lines identified in the FUSE spectrum. 
"*" denotes lines with possibly inaccurately known wave- 
lengths. ":" denotes uncertain identifications. 
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Consequently, the majority of the detected absorption lines 
were never observed before in any astrophysical object. 

3. Atomic data 

We have performed a detailed search for reliable atomic data 
in a number of atomic data sources. These are the following 
databases: 

(i) Kurucz (1991) level and line lists 1 , 

(ii) Opacity Project (OP, Seaton et al. 1994) TOPbase 2 , 

(iii) National Institute of Standards and Technology (NIST) 3 , 

(iv) Chianti database (Young et al. 2003) 4 , 

(v) Kelly Atomic Line Database 5 , 

(vi) University of Kentucky Atomic Line List 6 , 

as well as the monograph of Bashkin & Stoner (1975). The 
Kurucz and OP databases are our basic sources for atomic data 
of heavy and light metals, respectively. Bashkin & Stoner is an 
important source for accurate level energies and line positions 
for light metals. Important additional data in this respect were 
found in the NIST and Chianti databases. 

1 http://cfa-www.harvard.edu/amdata/ampdata/kurucz23/sekur.html 

2 http://legacy.gsfc.nasa.gov/topbase/home.html 

3 http://physics.nist.gov/ 

4 http://wwwsolar.nrl.navy.mil/chianti.html 

5 http://cfa-www.harvard.edu/amdata/ampdata/kelly/kelly.html 

6 http://www.pa. uky.edu/~ peter/atomic/ 



Table 3. Photospheric oxygen lines identified in the Chandra 
spectrum. ":" denotes uncertain identifications. OV lines are 
not detectable ("-") but included here, because their absence is 
discussed in the text. 
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For the light metals (C-Na) the most serious problem in 
many cases is the lack of accurate line positions for highly ion- 
ized species, often making individual line identifications dif- 
ficult or impossible. The same problem is faced with Ca and 
the iron group elements: Their combined opacity is important 
for the overall flux distribution, but identifying these species by 
means of individual lines is hardly possible. However, it turns 
out that the cumbersome work for a critical data compilation 
can yield satisfying results. A large number of multiplets from 
Ne and Mg can be identified in the Chandra data. 

Line broadening is a critical issue and introduces consid- 
erable uncertainties for detailed line-profile fits. Linear Stark 
broadening is assumed for ions with one valence electron (He n, 
C iv, O vi) and treated in an approximate way, while quadratic 
Stark broadening is assumed for all other ions (for details see 
Werner etal. 1991). 

Another problem is posed by the requirement of good pho- 
toionization (b-f) cross-sections for metals. The strengths of 
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Fig. 2. Depth dependence of temperature, electron density, and ionization fractions of chemical species in the model with 
r eff =175 000K. 



absorption edges within the Chandra spectral range and the ad- 
jacent EUV/FUV region affects the spectral energy distribu- 
tion. We have investigated this problem, as well as the effect 
of pressure ionization on the absorption edges, in detail before 
and refer to our earlier EUVE analysis of H 1504+65 (WW99). 
What is a new problem, that becomes important with the supe- 
rior Chandra spectral resolution, is the occurrence of autoion- 
ization resonances in the b-f cross-sections. Some of these reso- 



nances are strong and narrow, mimicking absorption lines in the 
spectra (Fig.^. Although these resonances are included in the 
OP data, their wavelength position is too uncertain for a unique 
identification in the observations. As a consequence, hitherto 
unidentified lines in the Chandra spectrum may be due to b-f 
resonances of metals and, vice versa, absorption features in 
the synthetic spectra without an observed counterpart could be 
due to such resonances in the models. For this reason we have 
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Table 4. Photospheric neon lines identified in the Chandra 
spectrum. ":" denotes uncertain identifications. 



A/A 




Transition 


Ne vi 


111.10, .16, .26 


2p 2 P° - 3p 2 D 


114.07, .13, .24, .30 


2p 2 P° - 3p 2 P 


120.21, .24, .27, .31, .35, .40, .45 


2p 24 P -3d 4 P° 


121.05, .06, 


.10, .11, .13, .15, .19, .21 
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2p 23 P - 3d 3 P° 


110.53, .56, 
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2p 23 P -3d 3 D° 


111.15 




2p 21 D - 3d'P° 


111.81 




2p 2 'D - 3d 'F° 


115.33, .39, 


.52 


2p 3 P° - 3s 3 S 


115.96 




2p 2 'D - 3d'D° 


116.69 




2p 'P° - 3d'D 


120.20, .27, 


.33, .35, .42, .48 


2p 23 P ~3s 3 P° 


121.13 




2p 2l D -3s'P° 


121.77 




2p 21 S -Sd'P" 


127.67 




2p 'P° - 3s 'S : 


133.64 




2p 21 S ~3s'P 


135.29, .33, 


.38, .40, .50, .54 


2p 23 P -3p 3 P° 


141.22 




2p 2 'D - 3p 1 P° : 


Ne vm 


88.08, .12 




2s 2 S - 3p 2 P° 


98.11, .26 




2p 2 P° - 3d 2 D 


102.91, 103.08 


2p 2 P° - 3s 2 S 



Table 5. Photospheric magnesium lines identified in the 
Chandra spectrum. ":" denotes uncertain identifications, "bi" 
denotes blends with lines of other elements. 



A/A 
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Mg v 


137.40, .41, .42, .74, .75, .8£ 
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qa oa 
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98.98, .99, 99.02, .04 


2p 3 2po 


- 3d' 2 P 




100.70, .90 


2p 3 2 D° 


- 3d 2 F 




101.49, .55 


2p 3 2 D° 


- 3d 2 P 


Ovi bl 


104.52, .53, .59, .60 


2p 3 2 P° 


-3d 2 P 


Ovibl 


111.17, .19 


2p 3 2 P° 


- 3s" 2 S 


Ne vn bl 


111.55, 75, .86 


2p 3 4 S° 


- 3s 4 P 




113.19 


2p 3 2 D° 


- 3s' 2 D 




116.97, 117.22 


2p 3 2 D° 


-3s 2 P 




116.97, .99 


2p 3 2po 


- 3s' 2 D 




117.55 


2p 42 D 


- 3s v 2 P° 




121.01, .03, .29, .30 


2p 3 2 P° 


- 3s 2 P 


Ne vn bl 


123.59 


2p 42 D 


- 3s iv2 D° : 




130.31, .64 


2p 42 P 


- 3s v 2 P° 




137.81, 138.17 


2p 42 P 


- 3s lv 2 D° : 





substituted the OP b-f data by hydrogen-like cross-sections in 
the detailed final spectral fits shown in this paper. This does 
change the continuum flux slightly, but in view of the difficul- 
ties we encounter with fitting the overall flux distribution (see 
below), this is considered unimportant. Generally, strong reso- 
nances are fortunately not very frequent and can be identified 
in the models on a case-by-case basis. 

Let us now present the design of our model atoms. A 
good orientation concerning the ionization degree of differ- 
ent species is provided by Fig.|2] Besides the temperature and 
electron density stratification of a model with 7^=175 000 K 
it shows the ionization fractions of all considered chemical 
species. We will discuss the changes in ionization balance with 
T eS in Sect.|5] 

3.1. Helium, carbon, and oxygen 

Because of the high temperatures in the atmosphere of 
H 1504+65, helium is strongly ionized so that a one-level Hei 
model atom is sufficient. 

For C and O, the helium-like noble-gas configurations (C v, 
Ovu) represent the highest ionization stages in the model 
atoms. They are also the dominant ionization stages in our 
model atmospheres (see Fig.|2ji. Test calculations based upon 
extended model atoms including H-like ions C vi, O vm showed 
that this is sufficient. The population of the latter ions be- 
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Table 5. continued 



A/A 




Transition 


Mg vii 



78.34, .41, .52 


2p 33 P 


-3p 


3po 




80.95, 81.02, .14 


2p 3 3 P 


-3p 


3 S o 




82.94, 82.97, 83.01 


2p 3 5 S° 


-3d 


5 P 




83.51, .56, .59, .64, .71, .76 


2p 33 P 


-3d 


3po 




83.91, .96, .99, 84.02, .09, .11 


2p 33 P 


- 3d 


3 D o 




85.34 


2p 21 D 


- 3d 


1 po 




85.41 


2p 2 'D 


- 3d 


Ipo 




87.72 


2p 2 'D 


- 3d 


'D° 




87.89 


2p 2 'D 


- 3d 


3po 


intercomb. 


88.68 


2p 21 S 


-3d 


1 po 




94.04, .17, .24 


2p 3 5 S° 


-3s 


S P 


Mg vm bl 


95.03, .04 


2p 3 3 D° 
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comes significant only in the outermost and innermost layers 
of the hottest models calculated (T e ff=250000K). This neither 
changes the observed C iv and O vi line profiles nor the flux dis- 
tribution in the wavelength bands studied here. Also, the pop- 
ulation of excited levels in C v and O vn is very weak, which 
explains why no subordinate lines from these ions are detected 
in the observed data (the respective resonance lines are located 
blueward of the recorded Chandra spectrum). 

The UV/optical spectrum of H 1504+65 is dominated by 
C iv and O vi lines (but note that no C iv lines are located in the 
Chandra spectral range). Available level energies often yield 
line positions that do not coincide with line positions from lab- 
oratory measurements. This is frequently the case for O vi lines, 



but also sometimes for C iv lines. When an exact line position 
could be taken from the databases, we shifted the synthetic line 
to that position, but uncertainties remain for many lines, par- 
ticularly in the FUV range. Similar problems occur with many 
O v lines in the Chandra range, preventing their identification, 
which is unfortunate because these lines would be sensitive T e g 
indicators. In all cases we restricted detailed profile analyses to 
lines with accurately known wavelengths. Tables [2] and [3] list 
the identified C and O lines. 



3.2. Neon, magnesium, aluminum, and sodium 

Besides O vi lines, lines from Ne vn are the strongest absorption 
features in the Chandra spectrum of H 1504+65. This ion also 
has a strong optical line at 3644A, which we detected in a Keck 
spectrum (WW99), and in the present paper we announce the 
detection of the Nevn 973. 3A absorption line 2p'P°-2p 2 'D; 
to our knowledge it is the first time ever seen in an astrophys- 
ical spectrum. There is some uncertainty in the literature con- 
cerning the accurate wavelength position of this line. The first 
measurement was published by Johnston & Kunze (1971) who 
give 973. 6A. Lang (1983) quotes 973. 33A (without reference) 
and this value is also found in the Chianti database. According 
to Kramida (NIST, priv. comm.) the best measurement for this 
line was done by Lindeberg (1972). The measured wavelength 
is 973. 302+0. 005A. For our synthetic spectra we have adopted 
973. 33A. We also model the line spectrum of Ne vi and Ne vm 
and successfully identify these ions in the Chandra data (see 
below). Table0]gives a complete list of all lines considered in 
the spectrum synthesis calculations. As already mentioned, fine 
structure splitting of multiplets is problematic in those cases, 
in which knowledge of energy levels within one term is incom- 
plete. This can cause uncertainties in wavelength positions of 
individual lines within a multiplet of the order 0.1 A. We have 
omitted all lines involving energy levels uncertain to an extent 
that the corresponding wavelength uncertainty is much larger 
than 0.1-0.2A. As a consequence, many unidentified lines in 
the Chandra spectrum probably stem from such Ne lines. 

After first test calculations with a small magnesium model 
atom we found that some lines of these species might be present 
in the Chandra data. However, because the Mg lines are much 
weaker than O and Ne lines, we have designed a very detailed 
model atom for Mg in order to look for as many lines of this 
species as possible. After a critical data compilation and the re- 
jection of many lines because of wavelength uncertainties there 
is no doubt that we can identify many lines of Mg v-vm. As in 
the case of Ne, many unidentified Chandra spectral lines cer- 
tainly stem from Mg lines with uncertain wavelength positions. 
It is interesting to note that we predict an observable intercom- 
bination line, which indeed is identified in the Chandra data 
(Table|5}. 

In the very same way we tried to identify Al and Na lines in 
the Chandra spectrum. Although large model atoms were con- 
structed, we cannot detect these species beyond doubt. Details 
of this search will be described below. 
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Fig. 3. The FUSE spectrum of H 1504+65 with the strongest photospheric lines identified. They stem from C iv and O vi. Note 
particularly the strong Nevn line at 973. 3A. This line is identified for the very first time in a stellar spectrum. All emission 
features are geocoronal, not photospheric. Overplotted is a model spectrum (r e ff=200000K). Observed and model spectra are 
convolved with Gaussians with FWHM=0.05A and 0. 1 A, respectively. 



3.3. Calcium and iron group elements 

Synthetic spectra including lines from Ca and Fe group ele- 
ments were calculated using a superlevel/superline approach 
which is described, e.g., in Rauch & Deetjen (2003). All line 
data were taken from Kurucz's line lists, except for Fe x, which 
were taken from the Opacity Project's TOPbase. Nix lines ex- 
ist neither in those two datasets nor in the NIST database and, 
hence, could not be considered here. Kurucz's data compila- 
tion consists of two sets. The first one includes only lines with 
accurate wavelengths determined from laboratory experiments 
("POS" data). The second list additionally includes predicted 
lines with uncertain wavelengths ("LIN" data). Our strategy is 
that we perform model and line formation iterations using the 
full (LIN) line list, because precise wavelength positions are 
unimportant at this stage of calculations, but accounting for the 
complete line opacities is important. The final synthetic spectra 
that are compared with observations, however, are computed 
using the POS data only, in order to facilitate line identifica- 
tions. We shall see that strong blanketing by the large number 
of LIN lines makes the identification of individual POS lines 
of iron and nickel virtually impossible. The OP Fe x lines are 
treated like Kurucz's LIN lines, because their positions are un- 
certain, too. 



4. Overview of FUSE and Chandra Observations 

H 1504+65 was observed by FUSE on four separate occasions, 
with a total integration time of 67 ksec. The FUSE instrument 
consists of four coaligned telescopes, each with a prime-focus 
spectrograph. The observation log is given in Table |6j all data 
were obtained in time-tag mode and processed with CALFUSE 
v2.4. The LWRS spectrograph apertures were used for all ob- 
servations, hence the zero-point of the wavelength scale is un- 
certain to within about +0.1 5 A. All exposures were photomet- 
ric, or nearly so, in all channels. Descriptions of the FUSE in- 
strument, and channel alignment and wavelength calibration is- 
sues are given by Moos etal. (2000) and Sahnow etal. (2000). 
The wavelength scales for spectra from the individual expo- 
sures were coaligned separately for each channel by means 
of narrow interstellar absorption lines. Spectra from channels 
other than LiFl were then given an additional zero-point offset 
so that H2 lines (which are present in all channels) had a com- 
mon heliocentric velocity. The spectra from all the channels 
were then resampled onto a common wavelength scale with 
0.025A pixels, and combined. The resulting spectrum covers 
912 - 1 187A with a typical spectral resolution of 0.08-0.10A. 

The spectrum is characterized by a few broad photospheric 
lines of C IV and O VI and many narrow interstellar lines, pre- 
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Fig. 4. Details from the FUSE spectrum. Left panel: region around the newly identified Ne vn line. Right panel: Close-up of the 
region near the He n (n=2— >7) line. It is not detectable. Overplotted are four models with increasing He abundance, confirming the 
He-deficiency of H 1504+65 (T e ff=200000K). Observed and model spectra are convolved with Gaussians with FWHM=0.05A 
and 0.1 A, respectively. 



Table 6. Log of FUSE observations of H 1504+65. 



Prog ID 


Date 


T cxp (ksec) 


Q 1090202 


2000-02-19 


11.3 


S60 10202 


2002-01-28 


32.7 


M1052601 


2003-01-13 


13.9 


M1052602 


2003-05-06 


8.9 



dominantly from H2 (see Hig-EJ. The identified photospheric 
lines are listed in Table |3 All other absorption features are of 
interstellar origin, except for two broad and shallow features at 
1163A and 1165A (marked by "?" in Fig.[3), which we think 
are of photospheric origin, but cannot be identified. They are 
also seen, even more prominently, in other very hot PG1159 
stars, the central stars Kl-16 and RXJ21 17+3412. There are 
S vm lines listed at these wavelengths, however, this identifica- 
tion is regarded as unlikely, because of the high level energies 
and because other lines from this ion should be detectable, too, 
which is not the case. The model spectra shown here are shifted 
by +0.09 A in order to fit the bulk of the observed lines. 

Two interesting details are shown in Fig.@] The left panel 
focuses on the region around a rather strong absorption line, 
which stems from Ne vn. As already mentioned this is the first 
detection of this line in a stellar spectrum. It is also seen in 
FUSE spectra of other PG1159 stars (e.g. in the prototype 
PG1 159-035 itself). The overplotted model, which fits the ob- 
served line very well, has a Ne abundance of 2%, in agree- 
ment with the results from the Chandra spectrum analysis de- 
scribed below. The right panel of Fig.@]is centered around the 
location of a He 11 line. He 11 absorption is not detectable and 
a comparison with models shows that H 1504+65 is strongly 
helium-deficient. This corroborates the result of an optical 
spectral analysis (He< 1 %), which was based on the lack of a 



He II 4686A emission line (W91). From Fig.|4]we derive a (less 
tight) upper limit for the He abundance, namely He<5%. 

H 1504+65 was observed with Chandra on Sep. 27, 2000, 
with an integration time of 7 hours. Flux was detected in the 
range 60A-160A and the spectral resolution is about 0.1 A. 
Fig.|5]shows the overall spectrum. It is characterized by a roll- 
off at long wavelengths due to ISM absorption. The maxi- 
mum flux is detected near 1 10A. Between 105 A and 100 A the 
flux drops because of photospheric absorption from the Ovi 
edge caused by the first excited atomic level. The edge is not 
sharp because of a converging line series and pressure ioniza- 
tion (see WW99 for detailed model spectra). Below 100 A the 
flux decreases, representing the Wien tail of the photospheric 
flux distribution. Fig.[5]demonstrates that our models cannot fit 
the entire wavelength range at a unique temperature, however, 
from the overall flux shape T e s can be constrained between 
175 000 K and 200 000 K. 

We have looked for possible explanations for the inability 
of our models to fit the overall X-ray flux distribution. We must 
look for processes which decrease the flux level at A > 100 A 
and/or increase the flux level at shorter wavelengths. Line blan- 
keting by heavy metals (introduced below) does not solve the 
problem. Including the large number of absorption lines blocks 
flux at A > 100A, but this flux re-emerges in the same wave- 
length region just in between spectral windows of strong line 
blanketing, hence, a flux increase at A < 100 A cannot be en- 
forced by including this opacity source. From this result we 
feel that we must look for a missing continuous opacity. 

One possibility was that our carbon model atom needed 
to be extended to higher ionization stages. Excited levels of 
Cv (2s and 2p) potentially cause absorption edges between 
130-150A. This could decrease the flux for A > 100A and at 
the same time, because of total flux constancy, could increase 
the flux A < 100A. However, a test calculation has shown 
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that this does not solve the problem. Even in a model with 
r e ff=250000K, in which C v is the dominant ionization stage 
of carbon in most parts of the atmosphere, only weak C v edges 
appear, decreasing the maximum model flux near 1 10A by less 
than 10%. By analogy, excited O vn levels can cause absorption 
edges. But their strength is expected to be even less important, 
because the ionization potential of O vi is higher than that of 
Civ. 

We also tried to vary other model parameters to overcome 
this problem. Increasing the surface gravity from log g — 8 to 
log g — 9 causes a flux decrease of about 50% in the 100-130A 
range. The additional opacity is caused by Ne vn and O vi line 
wings which become very strongly pressure broadened. In ad- 
dition, increased pressure ionization of O vi causes a shift of the 
strong ionization edges towards lower wavelengths. But again, 
the result is not a flux increase at A < 100 A, but a higher flux 
at A > 130A. A similar result was obtained from an exercise, 
where we increased the (uncertain) line broadening parameters 
even by very large amounts. The oxygen abundance is also not 
a critical parameter. Tests show that the strength of the O vi ab- 
sorption edge near 100 A is insensitive to variations of the O 
abundance within reasonable limits. Even a reduction of the O 
abundance by a factor of 10 decreases the continuum flux jump 
at the absorption edge from 95% only to 80%. In the same di- 
rection, we artificially decreased the b-f cross-section of the 
Ovi levels by factors larger than we expect from the uncer- 
tainty in the atomic data for this relatively simple one-electron 
ion. Switching between OP cross-sections and hydrogen-like 
cross-sections has only a weak effect on the model flux. 

To summarize these tests, the only way to significantly in- 
crease the short wavelength flux is to increase the effective 
temperature to roughly 200 000 K. An unknown opacity source 
must be responsible for suppressing the excess model flux at 
longer wavelengths. A detailed analysis of the Chandra line 
spectrum can be used to further constrain the effective temper- 
ature by using ionization balances of several species. We will 
now show that this supports the r e ff=200 000 K estimate. 

5. Analysis of the X-ray absorption line spectrum 

Figure |6] shows in more detail the Chandra spectrum of 
H 1504+65. Considering the problems with fitting the over- 
all continuum flux, we now normalize the model flux to the 
local continuum in each panel of this figure. In order to fa- 
cilitate line identifications, we overplot that model from our 
grid which turned out to fit best to the line spectrum. The ob- 
served spectrum was shifted by -0.06A to the rest wavelength. 
We detect lines from O vi, Ne vi-vm, and Mg v-vm. The strong 
and broad O vi and Ne vn lines were already detected in the 
EUVE spectrum (WW99), the detection of magnesium is new. 
Identification of even the strongest lines becomes difficult at 
longer wavelengths (see bottom panel of Fig. [6} due to the de- 
creased S/N-ratio. 

For a concise presentation of our analysis, we first concen- 
trate on the effective temperature which is the most important 
parameter determining ionization balances, while keeping fixed 
log g — 8 from our analysis of optical spectra. We think that an 
improvement of the error range for the gravity (+0.5 dex, W91) 



is not possible because of uncertainties in line broadening cal- 
culations. We also keep fixed the relative abundances of C, O, 
and Ne, which are the dominant species, and include He as a 
trace element (which has unimportant effects on model struc- 
tures and fluxes). Hence, the composition of all models shown 
here is: He=l%, C=48.5%, 0=48.5%, Ne=2%. As mentioned 
above, other elements are included later in line formation iter- 
ations. 

We start the T e ff analysis with a first estimate for the Mg 
abundance (2%) which will be justified later. After this we turn 
to the search for other species (Al, Na, Fe, Ni). 

5.1. Effective temperature from ionization balances 

Given the described uncertainty in the determination of J e ff 
from the overall X-ray flux shape, it is important to constrain 
the temperature by comparing line strengths from different ion- 
ization stages of the identified species. We discuss their ap- 
pearance in the context of three models with r e ff=175 000 K, 
200 000 K, and 250 000 K. Qualitatively, the shift of ionization 
balances and, hence, the change of line strengths with increas- 
ing r e ff can be understood with the help of Fig.|3 

Oxygen First of all, the observed spectrum does not exhibit 
Ov lines. Two relatively unblended Ov lines (at 119.1 A and 
124. 6 A) appear very strong in the 175 000 K model and de- 
crease in strength with increasing r e ff. Even at 200 000 K the 
lines are still quite strong. On the other hand, at 250 000 K all 
strong O vi lines generally become too weak. This clearly sug- 
gests r e ffw200000K). As we will see later, strong line blanket- 
ing by Fe group elements can severely hamper the detection 
of individual weak lines. Fig.[9] shows for example how the 
Ov 119.1 A line becomes difficult to detect in the respective 
model spectra. 

Neon We detect several Ne vi lines, the strongest are located at 
122.5, 122.7A and 1 14.1-1 14.4A. They come out too strong in 
the 175 000 K model and too weak at 250 000 K, they fit best at 
T eS =200 000 K. T eB =250 000 K can also be excluded because 
all Ne vn lines become too weak. 

Magnesium Increasing T eS in our models from 175 000 K 
to 250 000 K has the following consequences for the line 
strengths: Lines from Mg v and Mgvi become weaker, Mg vn 
strengths are hardly affected (they are maximal at 200 000 K), 
and Mgvm lines become very strong. We point out some 
individual lines, which altogether suggest 7^=200 000 K: 
Mgv 137.4A is barely detectable. It would have completely 
disappeared if 7^=250 000 K. A broad trough of Mg vn lines 
is located near 83-84A. It is very nicely fitted by our models, 
however, it is quite insensitive against T e ff. More interesting is 
Mgvm 74.3-74.4A. It is best fitted by the 200 000 K model, 
while too weak or too strong at 175 000 K and 250 000 K, re- 
spectively. It seems as if two lines from Mg vi are predicted too 
strong even at 200 000 K, namely the multiplets near 123. 6 A 
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Fig. 7. Synthetic sodium line profiles of three different ioniza- 
tion stages (one Navm and one Navi singlet and two Navn 
doublets) compared to the Chandra spectrum. There is no con- 
vincing evidence for Na lines in H 1504+65 (r eff =200000K, 
Na abundance 0.3% and 0.03%). Observed and model spectra 
are convolved with Gaussians with FWHM=0.02A and 0.03 A, 
respectively. Tick marks are spaced by 0.1 A. 

and 111 .7A but, again, Fe group line blanketing could obscure 
these features. 

As to the magnesium abundance, we calculated a 200 000 K 
model with the Mg reduced by a factor of 10. The vast majority 
of the lines become significantly too weak. In contrast the as- 
sumed abundance of 2% generally gives a good fit. We estimate 
the Mg abundance error to 0.5 dex. In view of the uncertainty 
to define the flux continuum we have to accept that a more ac- 
curate determination is not possible. 

5.2. Search for aluminum and sodium 

As already mentioned, the detection of weak individual lines 
is very difficult due to strong line blending from Fe group el- 
ements. Even careful modeling of this dense line curtain does 
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Fig. 8. Synthetic aluminum line profiles of two Alvn quar- 
tets (left two panels). The right panel shows two Alvi lines, a 
singlet (left line) and one triplet component. There is no con- 
vincing evidence for Al lines in H 1504+65 (r eff =200000K, 
Al abundance 0.1% and 0.01%). Observed and model spectra 
are convolved with Gaussians with FWHM=0.02A and 0.03 A, 
respectively. Tick marks are spaced by 0. 1 A. 

not help, because the positions of most Fe group lines are in- 
accurate. So we are very conservative in our search for further 
species in the Chandra spectrum. 

Sodium and aluminum are important elements for the dis- 
cussion of the evolutionary status of H 1504+65. Fig.[7]com- 
pares calculated sodium line profiles with the Chandra spec- 
trum. We have selected some of the strongest lines of three 
ionization stages, Na vi-vm. The profiles are computed for two 
different Na abundances, 0.3% and 0.03%. There is no clear ev- 
idence for sodium lines in the Chandra spectrum. We estimate 
an upper limit of Na=0.1%, which is about 30 times solar. 

Al vi and Al vn show the strongest aluminum lines in the 
model with T s g -200 000 K. Al v and Al vm are predicted to be 
much weaker. Some of these strongest lines are shown in Fig. [8] 
Two models were calculated with abundances of 0.1% and 
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Fig. 5. Overview of the Chandra spectrum of H 1504+65 (solid line). Two models with different T e g are shown. The 175 000 K 
model (long dashed line) fits the overall flux at long wavelengths, while it underestimates the short wavelength flux. The hotter 
model (200 000 K, short dashed line) fits at short wavelengths but overestimates the long wavelength flux. The model spectra 
were attenuated by an ISM model with «h = 9.02 • 10 19 cirT 2 and 1.06 • 10 20 cirT 2 , respectively (assuming «He/«H = 0.1, and 
"Hei/»Hen = 0.5), then folded through the instrument response and normalized to the observation to fit near 130A. The models 
include He, C, O, Ne, Mg as described in the text. For clarity, all spectra are smoothed with Gaussians (0.5A FWHM). 
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Fig. 6. Expanded view of the Chandra spectrum and a model with T e g =200 000 K (including He, C, O, Ne, and Mg as described 
in the text). 
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Fig. 9. Effect of line blanketing on the X-ray spectrum of H 1504+65. Top panel: only O and Ne lines included in the model. 
Middle panel: iron lines added. Bottom panel: all lines from elements Ca to Ni added. In the bottom panel all lines with calculated 
energy levels, i.e., uncertain wavelength positions, are included. This gives an impression of the strong line blanketing. It makes 
individual line identifications of heavy metal lines virtually impossible (7^=200 000 K). Observed count spectrum and model 
spectra (relative flux) are convolved with Gaussians with FWHM=0.02A and 0.03 A, respectively. 



0.01%. Although there seems to be a good fit by the calculated 
quartet at 75. 2 A (left panel), this is regarded as a chance co- 
incidence with some other unidentified feature, because other 
Al lines cannot be identified. We estimate an upper limit of 
A1=0. 1 %, which is about 20 times solar. 

5.3. Search for iron and nickel 

Figure [9] shows the strong influence of heavy metal line blan- 
keting. It becomes clear that the identification of individual 
lines from the POS line lists is affected by the many more lines 
from the LIN lists, which have uncertain wavelengths. Indeed, 
it is not possible to identify the majority of the strong Fe ix and 
Niix lines appearing in the synthetic spectra shown in Fig. 1101 
hence, we cannot derive abundances for individual Fe group el- 
ements. For all models including Ca-Ni we have assumed solar 
abundances for these elements. This appears to be a good esti- 
mate, because in many parts of the Chandra spectrum the line 
blanketed models seem to give a better overall fit to the obser- 
vation. So we conclude that the Fe group elements have a total 
abundance which is in accordance with solar values, albeit with 
a large uncertainty. 

5.4. Unidentified absorption features 

Many absorption features in the Chandra spectrum cannot be 
identified. Considering the large number of potential iron group 



lines whose exact wavelength position is unknown, this is not 
surprising. However, let us comment on about a dozen of the 
strongest of these features. Some of them can be attributed 
to magnesium and neon lines which are not included in our 
model atoms, because their high upper level energies would re- 
quire excessively large model atoms. These are lines from Mg v 
(1 14.8/1 15.1 A), Mg vii (79.25 A, 84.64A, 87. 13/. 17 A, 91 .49 A), 
Mg vm (80.89A, 84.92A), and Ne vi (109.03/.07A). The strong 
absorption at 125. 1 A remains unidentified. It cannot be a quar- 
tet transition of Mg vi located there. It is included in our mod- 
els but it is very weak. The same applies to a Nevi dou- 
blet at this position. A line observed at 86.4A could be from 
Alvm 86.43 A, but this is a triplet-quintet intercombination 
transition and therefore probably too weak to explain the obser- 
vation. Other strong unidentified features are located at 93. 5A 
(a Fe vm line here is too weak) and 1 13.4A. We recall that they 
could also stem from resonances in bound-free cross-sections. 

It is our general impression, supported for example by 
Fig. [9] that, besides strong O, Ne, and Mg lines, the blanket- 
ing by many relatively weak iron group lines characterizes 
the Chandra spectrum. Our models show similar characteris- 
tics, except for wavelengths below, say, 100A. We attribute this 
to those high ionization stages, particularly of Fe and Ni, for 
which no line lists are available. We are probably neglecting 
a large number of Fexi and Nix-xi lines at A < 100A, hence, 
atomic data for these species are badly needed. In the case of 
Ni x we computed the positions of a dozen of resonance lines 
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from Bashkin & Stoner (1975) energy levels. They are located 
around 85A and one of these lines could in fact explain the 
unidentified 86.4A feature mentioned above. 

6. Summary of spectral analysis 

Let us summarize the properties of H 1504+65: 

r eff = 200 000 K+ 20 000 K 
\ogg = 8.0 ±0.5 [cgs] 

Element abundances in % mass fraction: 

C = 48 
O = 48 
Ne = 2 
Mg = 2 
Fe - group = 0.14 (solar) 
He < 1 
Na < 0.1 
Al < 0.1 

The value of log g and abundances for C, O, and Ne were taken 
from previous work (W91, WW99). Estimated errors for abun- 
dances are: +20% for mass fraction of C and O, and a factor of 
3 for Ne, Mg, and Fe-group. 

Stellar mass and luminosity can be derived by comparing 
the position of H 1504+65 in the g-T e g diagram with theoretical 
evolutionary tracks. We use the post-AGB tracks of Blocker 
(1995) and derive: 

M/M Q = 0.8363$ 
logL/L = 2.45!^ 
d/kpc = 0.67^3 

Note that the mass of H 1504+65 is considerably higher than 
the mean mass of the PG1159 stars (0.6 Mq). The spectro- 
scopic distance was obtained by comparing the measured vi- 
sual flux (V=16.24, Nousek et al. 1986) with the flux of the fi- 
nal model (r eff =200000K, logg = 8): FL [5400 A] =3. 42 • 10~ 3 
erg/cm 2 /s/Hz). Interstellar reddening was neglected for this de- 
termination, because it is very low. In fact, the best model fit 
to the continuum shape of the UV/FUV spectrum taken with 
HUT provided E(B-V)=0 (Kruk & Werner 1998). 

7. Discussion 

PG1159 stars are extremely hot post-AGB stars, many of 
them have already entered the hot end of the WD cooling se- 
quence. Generally their atmospheres are H-deficient and pri- 
marily composed of He, C, O, and Ne (typical abundances are 
33%, 50%, 15%, 2%; Werner 2001). It is now generally ac- 
cepted that they expose intershell matter, i.e., material that was 
located between the H- and He-burning shells in the former 
AGB precursor star. In contrast to "normal" H-rich post-AGB 
stars the PG1 159 stars have suffered a late He-shell flash, dur- 
ing the post-AGB or WD stage. The flash induces mixing of the 



intershell (which has a mass of about 0.0 1M Q ) with the much 
less massive H-rich envelope (about 10~ 4 M o ). Quantitatively, 
the intershell element abundances (He, C, O, Ne) in stellar 
evolutionary models are now in good agreement with obser- 
vations. The observed high O abundances require strong con- 
vective overshoot in order to dredge up oxygen from the stellar 
C/O core into the convective intershell (Herwig et al. 1999). 

The high amount of 22 Ne in the intershell is generated by 
ff-captures on 14 N from CNO cycling (see, e.g., Iben 1995, 
or Lattanzio 2003). Provided the temperature is high enough 
(which is the case in relatively massive AGB stars) another a- 
capture on 22 Ne yields 25 Mg or 26 Mg by (a,n) or (a,y) reac- 
tions, respectively. This can explain the high amount of Mg that 
we found in H 1504+65, which is the most massive PG1159 
star. We also would expect to find a Mg enrichment in other 
PG1159 stars, but the lack of appropriate spectral lines in the 
optical/UV prevents its detection. 

Besides the production of Ne and Mg in the intershell 
of AGB stars during the usual He-shell flashes (thermal 
pulses), one expects processing of these elements to Na and 
Al if they are dredged up into the H-burning shell, namely 
through the Ne-Na and Mg-Al cycles: 22 Ne(p,y) 23 Na and 
25 ' 26 Mg(p,y) 26 ' 27 Al. The high amounts of Ne and Mg in 
H 1504+65 on the one hand and the fact that we could not 
detect strong overabundances of Al and Na on the other hand 
suggests that these cycles were unimportant for the creation of 
today's photospheric chemistry of H 1504+65. 

It is unfortunate that we are not able to determine quan- 
titatively the abundances of individual Fe peak elements. We 
expect that we should see effects of s-processing on these ele- 
ments, shifting away their relative abundances from the solar 
values. In fact, it has been found that PG1159 stars are Fe- 
deficient, and we suspect this is indeed a result of n-captures 
on Fe seeds (Miksa et al. 2002, Werner et al. 2003b, Herwig 
etal. 2003). 13 C is the main neutron source and it is cre- 
ated by 12 C from He-burning and hydrogen from the surface 
layer: 12 C(p,y) 13 N(y6 + v) 13 C(a;,n) 16 0. Another source can be the 
above mentioned 22 Ne(a,n) 25 Mg reaction. 

But how can we explain the He-deficiency of H 1504+65? 
Is it just an extreme PG1 159 star in which He was completely 
burnt up to C and O during the late He shell flash? This can only 
be a speculation. There are no appropriate model calculations 
supporting such a scenario. Another idea which we already 
discussed before (W91, WW99) is that H 1504+65 has an evo- 
lutionary history completely different from the other PG1159 
stars. H 1504+65 could have gone through carbon burning, i.e., 
it would now have a O-Ne-Mg core and we could see the C/O- 
rich envelope on top of it. Hence, H 1504+65 might have been 
one of the "heavy-weight" intermediate-mass stars (8 MqSMJ 
10 Mq) which form white dwarfs with electron-degenerate O- 
Ne-Mg cores. Evolutionary models by Iben etal. (1997) pre- 
dict high Ne and Mg abundances for the C/O envelope, also 
in agreement with our analysis. So it is not possible to decide 
if we are looking onto a C/O core or a C/O envelope unless, 
in the latter case, we should see direct evidence for C-burning. 
One such piece of evidence would be a strong sodium ( 23 Na) 
enrichment. The 23 Na abundance at the bottom of the C/O en- 
velope is comparable to that of neon (main isotope 20 Ne) and 
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Fig. 10. Synthetic spectra (r ef j =200 000 K) with Fe ix lines (top panel) and Ni ix lines (bottom panel) compared to the Chandra 
spectrum. The POS line lists are used here to enable line identifications. Identifying individual lines is impossible, probably 
due to strong line blending by other heavy metal lines, as demonstrated in Fig.|5] Observed count spectrum and model spectra 
(relative flux) are convolved with Gaussians with FWHM=0.02A and 0.03 A, respectively. 



magnesium ( ' ' Mg, see Fig. 34 of Iben etal. 1997). This is 
not observed in H 1504+65, however, we stress the difficulty 
of giving a tight upper abundance limit for Na from the present 
observational data. 

The possibility that H 1504+65 is a O-Ne-Mg white dwarf 
is remarkable, because evidence for the existence of such ob- 
jects is rather scarce (Weidemann 2003). Evidence from single 
massive WDs is weak, and the most convincing cases are WDs 
in binary systems. Strong neon overabundances are found in 
novae (Livio & Truran 1994) or in eroded WD cores in low- 
mass X-ray binaries (e.g. Juett etal. 2001, Werner etal. 2004). 

It is also interesting to note that the supersoft X-ray source 
RX J0439.8-6809 was suggested to be similar to H 1504+65. 
This very faint object (V^2l) has a steep blue and almost fea- 
tureless optical/UV spectrum, except for each two Ovi and 
N v emission lines in optical spectra (van Teeseling et al. 1999, 
Reinsch et al. 2002). The presence of nitrogen in completely 
3o , -processed matter is hard to understand and, instead, it might 
be possible that the respective optical lines stem from C v and 
O vin. Such highly ionized emission lines have been detected in 
the hottest helium-rich DO white dwarf KPD0005+5106 and 
some PG1 159 stars (Werner et al. 1994) and were attributed to 
non-photospheric shock-heated emission regions. 

8. Summary 

We have analyzed new FUV and soft X-ray spectra of the 
unique object H 1504+65. We confirm its exotic chemical com- 
position, which is dominated by C and O. We confirm the 
high Ne abundance and find a similarly high abundance of Mg. 
This chemistry either reflects that of the core of a C/O white 
dwarf or the C/O envelope of a white dwarf with a O-Ne- 
Mg core. It therefore remains unclear if H 1504+65 has gone 
beyond 3a burning through a subsequent C burning phase or 



not. In any case, the origin of the He-deficiency remains ob- 
scure. H 1504+65 could be an extreme PG1159 star which - 
in contrast to the other stars of this group - has for some un- 
known reason burned up its helium completely. Alternatively, 
H 1504+65 could have burned carbon, now being a O-Ne-Mg 
white dwarf. Some unidentified mechanism (C shell flashes?) 
may be responsible for the loss of helium by ingestion and 
burning in deep hot layers. Interestingly, Iben and collabora- 
tors have predicted that such C -burning stars in the super- AGB 
phase could loose their H- and He-rich envelopes by a radia- 
tion driven superwind (Ritossa etal. 1996). H 1504+65 might 
resemble the result of such a scenario. 

A surprisingly rich photospheric absorption line spec- 
trum in the soft X-ray regime has been revealed by our 
Chandra observation. Although the overall flux distribution 
cannot be explained by a single model with a particular tem- 
perature, the ionization equilibria of O, Ne, and Mg suggest 
that r e jf is slightly higher than determined in previous analy- 
ses (200000K+20000K). This makes H 1504+65 the hottest 
known post- AGB star and white dwarf ever analyzed in detail 
with model atmosphere techniques. 
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